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Abstract 

We study the rare decays of B s — >■ M s il (M s = r],r]',(p and £ = u, e,/i, r) in the stan- 
dard model. With the hadronic form factors evaluated in the light front and constituent 
quark models, we estimate the decay branching ratios, respectively. We also discuss the 
longitudinal lepton polarization and forward-backward asymmetries. 



1 Introduction 



The inclusive flavor- changing neutral current (FCNC) process of B — > X s £ + £~ , which is 
suppressed and induced by electroweak penguin and box diagrams in the standard model 
(SM), has been recently observed by the Belle detector at the KEKB e + e~ asymmetric- 
energy collider [1] with the branching ratio (BR) of (6.1 ± 1.4lj;;f) x 1CT 6 for di-lepton 
masses greater than 0.2 GeV, where £ is either an electron or a muon and X s is a hadronic 
recoil system that contains a kaon. The exclusive decay of B — > K£ + £~ has been measured 
with the BR of (0.75±g;lf ± 0.09) x 10~ 6 [2] and (0.78lg;^18;iJ) x 10~ 6 [3, 4] at the Belle 
and Babar detectors by using 29.1 fb^ 1 and 77.8 fb~ x data samples, which agree with 
the theoretically estimated values [5, 6, 7, 8], respectively. Experimental searches at the 
B-factories [2, 3, 4] for B — > K*£ + £~ are also close to the theoretical predicted ranges 
[5, 6, 7, 8]. At Babar, an excess of events over background with estimated significance 
of 2.8 a has been observed and BR[B -> K*£ + £~) = (1.68lg:il ± 0-28) x 10~ 6 has been 
obtained [4]. It is clear that these FCNC rare decays are important for not only testing 
the SM but probing new physics. 

In this paper, we study another type of the exclusive decays, B s — > M s £ + £~ (M s = 
7], 7]', 0), related to the transition of b — > s£ + £~ at the quark level as shown in Figure 1. The 
decays of B s — > r]£ + £~ have been estimated in Ref. [9] and B s — * 4>£ + £~ have been explored 
in Refs. [10, 11]. Recently, the Collider Detector at Fermilab (CDF) Collaboration has 
given an upper limit of BR(B S — > 0/i + /i~) < 6.7 x 10~ 5 at a 95% confidence level [12]. In 
future, more data on the B s decays will be collected in various hadron colliders such as 
CDF-II, HERA-B and ATLAS. To study the decay rates and branching ratios, we need 
to calculate the transition form factors of the vector, axial-vector and tensor currents. 
There are many different candidates for this purpose, e.g., lattice QCD [13], QCD sum 
rule [14], and relativistic quark models. In this work, we use the framework of the light 
front quark model (LFQM) [15, 16] to evaluate the form factors and compared with those 
in the constituent quark model (CQM) [17]. 

This paper is organized as follows. In Sec. 2, we calculate the form factors for 
B s — > 77,77' and <fi transitions in the LFQM. We also present the results with those in 
the CQM. In Sec. 3, we study the decay rates and asymmetries of B s — > P(V)£ + £~ with 
I — v,e, /i,t and P(V) = pseudoscalar (vector) mesons. Our conclusions are given in Sec. 
4. 
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2 Form Factors 



2.1 Effective Hamiltonians and Matrix Elements 

The contributions to the decays of B s — > M s ££ (M s = rj, r/', 0) in the SM arise from the 
W-box and Z{pj)- penguin diagrams as seen in Figure 1. The effective Hamiltonians of 
b — > svv and b — > sl + l~ are given by [18] 

G p a - 

H = "7^ 9 B i^a XtD ( x *) hl » ~ 75 ^ s ^ ^ ~ 75 ^ ^ 
y 2 2tt sin (714/ 



and 
/X = 



C 8 e// (/i) SLl,b L hn + C 9 s l7 ^l £/W - 2mb ^^ s L ia^b R l-fl 



, (2) 



respectively, where Xt = m 2 /m 2 v , X t = V t *Vtb, D(x t ) is the top-quark loop function [19, 20] 
and Cs(fj), C$ and CV(/i) are Wilson coefficients (WCs) with their explicit expressions 
given in Ref. [20] for the SM. We note that C 9 is free of the /i scale. Besides the short- 
distance (SD) contributions, the main effect on the decays is from cc resonant states such 
as ^ and i.e., the long-distance (LD) contributions. To including the LD effect, in 
Eq. (2) we replace C 8 (/j) by C 8 e// (/i) [8, 20], given by 

C e 8 ff („) = CM + (3d (/,) + C 2 (/i)) (h (x, s) + I £ k Sj^J+n^ ) , (3) 

where we have neglected the small WCs, and /i(x, s) describes the one-loop matrix ele- 
ments of operators 0\ = s a ^i^P L bp c^^P L c a and 2 = s^^P^ c^^Plc [20], Mj (Tj) are 
the masses (widths) of intermediate states, and kj = — 1/ (3Ci (//) + C 2 (//)) [8]. The WCs 
at the scale of /x ~ m fe ~ 4.8 GeV are shown in Table 1. 

Table 1: Wilson coefficients for m t = 170 GeV and \x = 4.8 GeV. 



WC d 


c 2 


C 7 


c 8 


Cg 


-0.226 


1.096 


-0.305 


4.186 


-4.559 



To get the transition matrix elements of _B S — > M s with various quark models, we 
parametrize them in terms of the relevant form factors as follows: 



(77 S ( P2 )| V, \B s ( Pl )) 
(Vs(P2)\ T^ v q u \B s {p 1 )) 



F + (g 2 )P M + F_(g 2 )^. 



1 



m Bs + m r) ( r) / 



q 2 P,-(P-q)q, F T {q 2 ) , 



2 



(<f>{p2,e)\ IBM) 



(0(p 2 ,e)| (T fU ,±T%,) l f\B a (p i 



i - 



\-iV{q 2 )s, ua(5 6* v P a qP 

m Ba + L 

±A (q 2 )(P-q)e;±A + (q 2 )(e*- Pl )P, 
±A_(q 2 )(e* ■ Pl )q, 



±a (q 2 )(P-q) 
±a + (q 2 )(e*- Pl ) 



(4) 



where ttim = -B s , 77, 77' and 0) are the meson masses, ^1(^2) is the momentum of 
the initial (final) meson, e is the polarization vector of the vector meson <p, P = pi + P2, 
q = Pi- P2, = q 2 ^qi, A^ = 927^7551, T v» = ^icr^qi, T% = <Z2«<V75<7i, and F± ;T , V, 
^o,±5 9-, an d a o,± are the form factors. 

The daughter mesons in Eq. (4) are in the bound state ss with the physical masses, 
and rriej,. In the transfer matrix elements of B s — > r] s , there are two sets of the 
form factors, F± ;T , corresponding m v and m v /, respectively. 

The pseudoscalar meson 77 and i]' composed by the states mm, dd and ss can be ex- 
pressed as the mixtures of two orthogonal states 77^ and ?7 S , given by 77^ = (uu + dd)/y/2 
and 77s = ss [21]. Explicitly, one has that 



cos 9 — sin 9 
sin 9 cos 6* 



with 9 = (39.3° ± 1.0°) [21]. To study the decay rates of B s 
consider B s — > r\ s {m^)££ and _B S — > r) s (m v /)££ and then use 

r(5 s -> 77^) = sin 2 r[5 s -> 77 s (m^)££] 
T(fi s -> 77'^) = cos 2 r[B s -> 77 s (m T/ )^] . 

2.2 Light Front Quark Model 



(5) 

(rj,T]')££, we will first 



(6) 



Since the calculations of the transition form factors in Eq. (4) belong to the nonpertur- 
bative regime, the phenomenological quark models may be needed. One thing worthwhile 
mentioning here is that all of form factors will be studied in the time-like physical meson 
decay region of < q 2 < {m Bs — m VtV >^) 2 . As q 2 decreases, corresponding to the increasing 
recoil momentum, we have to start considering relativistic effects seriously. In particular, 
at the maximum recoil point of q 2 = where the final meson can be highly relativistic, 
there is no reason to expect that the non-relativistic quark model is still applicable. A 
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consistent treatment of the relativistic effects of the quark motion and spin in a bound 
state is a main issue of the relativistic quark model. 

The LFQM [22, 23] is the relativistic quark model in which a consistent and fully 
relativistic treatment of quark spins and the center-of-mass motion can be carried out 
[24]. The LFQM has been applied to study the heavy-to-heavy and heavy-to-light weak 
decay form factors in the timelike region [15, 25]. These calculations are based on the 
observation [26] that in the frame where the momentum transfer is purely longitudinal, 
i.e., q± = 0, q 2 = q + q~ covers the entire range of momentum transfers. The price one 
has to pay is that, besides the conventional valence-quark contribution, one must also 
consider the non-valence configuration (or the so-called Z graph) arising from the quark- 
pair creation in the vacuum. Unfortunately, a reliable way of estimating the Z graph is 
still lacking. However, the non-valence contribution vanishes if q + = 0, and it is supposed 
to be unimportant for heavy-to-heavy transitions [15]. In this paper, all of the values 
obtained from the LFQM are based on the formulas in Refs. [15, 16]. 

A meson bound state consisting of a heavy quark qi and an antiquark q 2 with total 
momentum P and spin S can be written as 

i"<™*» - /ww W(p - ft - fe) 

x £ * ss *(pi,p2,Ai,A 2 ) |gi(pi,Ai)g 2 (p 2 ,A 2 )), (7) 
where p\ and p 2 are the on-mass-shell light front momenta, 

P = (p + ,P±) , P± = (p\p 2 ) , P~ = — ~j (8) 

with 

p+ = (1 - X )P + , p+ = xP + , 

Pi± = (1 - x)P± + k±, p 2 ±=xP ± -k±. (9) 

Here (x, k±) are the light-front relative momentum variables, and k± is the component of 
the internal momentum k = (k±,k z ). The momentum-space wave-function ty SSz in Eq. 
(7) can be expressed as 

* ss '(pi,P2, Ai, A 2 ) = R s x %(x,k ± ) <j>(x,k ± ), (10) 

where (f>(x, k±) describes the momentum distribution of the constituents in the bound 
state 



with u being the meson scale parameter and R\f\ 2 constructs a state of definite spin 
(S,S Z ) out of light-front helicity (Ai,A 2 ) eigenstates. Explicitly, in practice it is more 
convenient to use the covariant form for R\f\ 2 [27]: 



R S x %(x, k ± ) = u( Pl , XjTvfa, A 2 ), (12) 



where 



J 75 (pseudoscalar S=0) 



and 

M = v/Mq 2 - ( mi - m 2 ) 2 (14) 

with 

2 _ m 2 + A; 2 m 2 + A; 2 
M ° " (1-x) + x (15) 

We normalize the meson state as 

<M(P', 5', S' g )\M(P, S, S z )) = 2{2-nfP + 5\P' - P)5 s >sSs>s z , (16) 

so that the normalization condition of the momentum distribution function can be ob- 
tained 

dx d 2 k i 



/ 



(*,*±)r = i- (it) 



2(2tt) 3 

We note that the form factors in Eq. (4) depend on the meson (M = (71(72) wave 
functions ^m^, To fix the parameters in the wave functions, one may use the meson 
decay constants /m, given by 

fu = V2A[^p^V M (x,k ± ) A (18) 
J 2(2tt) 3 ^ A 2 + k 2 ± 

where A = m qi x + m q2 (l — x) with m qi being the composed quark masses. 

In our numerical study, we use fs a = 200 MeV, and f Vs = (1.34 ±0.06) f n [21] which 
is the decay constant of the constituent bound state ss. By putting these decay constants 
into Eq. (18), we get the adapted u value for the wave function in Eq. (11). Explicitly, 
we find that uj Bs = 0.56 and uo ris = 0.45 GeV. 

In Figure 2, we display the form factors as functions of s = q 2 /M^ s for the transitions 
of B s — > i] s (with physical masses m Vtrj r) and B s — > <p with the parameters in Table 2. In 



Table 3, we show the form factors of the pseudoscalar daughter meson with q 2 = 0. Since 
F + does not depend on the daughter meson mass, we find that F + has the same value 
for both m v and m v /, whereas F_ and Ft are different. The form factors for B s — > <fi 
including vector and tensor current matrix elements are given in Table 4. We note that 
we have used the same scale parameter for the pseudoscalar and vetor mesons. In Tables 
3 and 4, we have also listed the results in the CQM [17] as comparisons. 

We note that, as shown in Figure 2, there are deviations for the form factors, in 
particular Ft, between the LFQM and CQM for the region of large q 2 values. As we 
know that in the LFQM the form factors are calculated in all physical q 2 region, some 
form factors, such as Ft, could increase as q 2 near the zero recoil point since the matrix 
elements depend on the overlap integral of the initial and final meson wave functions. 
While in the CQM [17], all form factors have similar q 2 behaviors since they have the 
same double pole forms of 



W) = I (19) 

Fi(0) l + a lS + a 2 s 2 v ; 

where s = q 2 /m 2 Bs , Fj(0) are the form factors at q 2 = 0, and <7 lj2 are the fitted parameters. 
Table 2: The parameters used in the calculation of the form factors (in GeV). 



m Bs m v m v i 




m h 


m s f B , 


fvs 


w Bs w Vs 


5.3696 0.5473 0.9578 


1.0194 


4.85 


0.35 0.200 


0.183 


0.56 0.45 



Table 3: Form factors F± t T for B s — > r) s (m ri ,m ri ') transitions at q 2 = in the LFQM and 
CQM. 





B s 


-> r] s (m v ) 




B s 


-> r]s(m, 


/) 




F+ 




Ft 


F + 




Ft 


LFQM 


0.354 - 


-0.360 - 


-0.369 


0.354 - 


-0.324 


-0.404 


CQM 


0.357 - 


-0.304 - 


-0.365 


0.357 - 


-0.304 


-0.390 
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Table 4: Form factors V, A 0j ±, g and a ,+ fc> r B s — > for P s — > transitions at g 2 = in 
the LFQM and CQM. 



P s ^0 




V 


A 


A + 


A_ 


9 


a 


a + 


LFQM 


-0.440 


-0.464 


0.276 


-0.295 


0.377 


0.374 


-0.374 


CQM 


-0.445 


-0.506 


0.310 


-0.379 


0.380 


0.380 


-0.380 



3 Decay Rates and Asymmetries 
3.1 Decay Rates 

Using Eqs. (1), (2), and (4), we obtain the decay amplitudes of B s — > (rj s , 4>)ll as follows: 



Gl °"" V ' D^F+P^py (1-75)1/, 



2V2 27rsin 2 W 



{ m Bs + 777.0 



2\/2 27rsin 2 W ~ 
+A (P-g)e; + A + ( e -P)Pj } i?y(l- 75 )„ , 



(20) 



(21) 



M(P s ^7^ + 



-M(P S 



GF^ern^t 



U 2m^T(g 2 ) |P|| 



£7^ 



x 



+ [C 9 F + P M + C 9 F_ ?M ] £VW} , 

G\FCH em A t 2 

2^^ m Bs 

1 G v e, uaP e*"P V + G° A (P- q) e; + G+ (e* • P) P, 
iiVw e w PV + F°(P.g)e; 



(22) 



£7^ 



+ 



where the functions Gy, Fy, G^, F A , G\ 



+ F+(e*.P)P /1 + F j 7(e*-P)< Z/1 

G\, Ft, and Fj are 



C 8 e// V 


C 7 m b g 


2(1 + ^) 


s 


C 8 e// A, 




2(1 + 


s 




C 7 rhba + 


2(1 + V^) 


s 


C 9 A_ 





F v = 



i defined by 

c 9 V 



(23) 



2(1 + 



C 9 A 
l + y/r, 
C 9 A + 



(24) 
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with s = q 2 /m 2 Bs , m b = m b /m Bs and = m 2 /m 2 Bs . From Eqs. (1), (2), (6) and (20)-(23), 
the differential decay rates for B s — > M s vv are found to be 



dT (B s -» M s vv) G 2 F \\ t \ 2 a 2 em \D{x t )\ 2 m\ 3 1/2 
¥m s J~m b 



where 



ds 

F v = sin 2 6 x ip v \F + \ 2 , 
T^i = cos 2 9 x {p v , \F + \ 2 , 



2 8 tt 5 sin 4 9 



'w 



•F. rh = 3s 



(l-v^,) 2 |Ao| 2 + 



(1 + y/?J 



IA)I 



4^ '" Ul (1 + v^) 2 
(l-r>-s)(l- y/^) 

+ 2r,(l + VF,) 

For B s -> M s £+£" (£ = e, //, r), we get 



|vf + 



4r>(l + v/F 



dr(B s - M s £+r) _ G||A t | 2 m| s a e 2 m 



3 • 2 9 vr 5 



/, 4t l 1/2 



, 2t\ 

1 + — I a Ms + t o 



Ah 



where 



t = mi /m B , r Ms = m M Jm B , 



Vm s = (1 - ?"AfJ 2 - 2s (1 + r Ms ) + s 2 . 
The formulas of and <5m s are given by 

= sin^x^f|^F + --^^| 2 + |C 9 F + | 2 ') , 



Oir, 



(25) 



(26) 



(27) 



(28) 



sin 2 9 x 6|C 9 | 2 {[2 (1 + r v ) - s] \F + \ 2 + 2 (1 - r v ) Re(F + F*) + s|F_| 2 } , (29) 



a„i = cos 9 x 



S v , = cos 2 ^x6|C 9 | 2 {[2(l + r v )-s]|F + | 2 + 2(l-r v ) J Re(F + F:) + s|F_| 2 } , (30) 



4s 



3(1-^) 2 (|G°| 2 + |^| 2 ) + 2^(|GH 2 + |F, 



1 2\ 



r <$> 



(i-^) 2 (|G°| 2 + |^| 2 ) + M|G+| 2 + :/;;ri 

+2(1 - r>)(l - r> - s)i?e(G'^G'+* + 
-48^|F y | 2 -72(l-r^) 2 |F°| 2 + 
12(1 -r^) 



■ a -72(l-r,)W+ 6[2(1 + ^" al ^ |a 
+ ?V^l| 2 + 



6s 
— < 

r 4> 



r <t> 



-^Re(F° A F+* + F° A FX* + F+F A 



(31) 



Table 5: Decay branching ratios of B s — > (77, r]', <p)uu. 







LFQM 


CQM 


10 6 Br( J B s - 




2.34 


2.17 


lO^rfs, - 


-> 77'1/P) 


2.52 


2.38 


K^Brf.B, - 




12.02 


11.65 



Table 6: Decay branching ratios of B s — > (rj, rf, <p)£ + £ without including LD effects. 

without LD 

Decay Mode LFQM CQM Ref. [9] Ref. [10] Ref. [11] 
10 7 Br( J B s r]e + e-) 3.43 3.13 
10 7 Br(^ s -> ?7// + /0 3.42 3.12 4.6,5.2 
10^(5, -> ^r+r-) 0.67 0.67 
10 7 Br(S 8 -> r/e + e~) 3.64 3.42 
10 7 Br(S s -> 77W") 3.63 3.41 
10 7 Br(5 s -> 77'r+r") 0.50 0.43 

10 6 Br(£ s -> 0e+e~) 1.72 1.69 - - 2.01,1.87 

10 6 Br(B s -> 1-64 1.61 - 2.5 1.65,1.25 

10 7 Br(5 s -> 0r+r-) 1.60 1.51 - - 1.38,2.28 



By using the form factors of the LFQM and CQM in Figure 2 and |A t | = |T4 6 T4 S | ~ 
0.041, from Eqs. (25) and (27), we now can estimate the numerical values for the decay 
rates. In Figures 3 and 4, we show the differential decay branching ratios as a function 
of s = q 2 /m 2 Ba for B s — > M s ££ (£ = u, e, //, r) in the two models, respectively. In Tables 5 
and 6, we summarize the BRs of B s — > (rj, rj', (fi)£ + £~ (/ = u, e, //, r), where LD effects for 
the charged lepton modes are not included. In Table 6, we also show the results in Refs. 
[9], [10] and [11]. In the table, the two values under Ref. [9] represent two inputs of the 
r) mixing parameter and those under Ref. [11] are from two kinds of hadronic models. 

With the LD effects, we need to introduce some cuts around the resonances of J ftp 
and ip'. Explicitly, we use the same cut as the one in the CDF experiment [12], which 
leads to three regions given by 

I. 2m e < yg2 < 2.9; 

II. 3.3 < yg^ < 3.6; 

///. 3.8 < < m Bs -m Ma . (32) 

In Table 7, we present the decay BRs in terms of the regions shown in Eq. (32). From 
Table 7, it is interesting to pointed out that the results of Br(B s — > (f)fi + fi~) = 1.23(1.24) x 
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1CT 6 in the LFQM (CQM) are about two orders of magnitude less than the recent CDF 
experimental upper limit of 6.7 x 1CT 5 [12]. We note that with the LD effects, in Refs. 
[10, 11] the branching ratio of B s — > is found to be around 1.9 and (1.91, 1.45) x 

10~ 6 , respectively. The values for Ref. [11] were derived with a much smaller cut. 

Table 7: Decay branching ratios of B s — > (rj, rj', 4>)£ + £~ with LD effects and the cuts. 



with LD 



I 



II 



III 



I+II+III 



Decay BR 



LFQM CQM 



LFQM CQM 



LFQM CQM 



LFQM CQM 



10 7 Br(£ s 

10 7 Br(B s 

10 8 Br(S s 

l0 7 Br(B s 

10 7 Br(5 s 

10 8 Br(£ s 

10 7 Br(5, 

10 7 Br(fi s 

10 8 Br(S s 



► r/e T e ) 

■ T)T + T~) 

■ r)'e + e~) 

in + n~) 

t)'t + t~) 

<l>n + n~) 

■ (f)T + T~) 



1.86 
1.85 

2.46 
2.44 

8.71 
7.91 



1.81 
1.80 

2.42 
2.40 

9.04 
8.30 



0.32 
0.32 
0.13 
0.34 
0.35 
0.14 
1.88 
1.88 
0.48 



0.28 
0.28 
0.12 
0.30 
0.31 
0.13 
1.83 
1.83 
0.48 



0.61 
0.61 
4.95 
0.27 
0.27 
3.27 
2.56 
2.56 
8.87 



0.49 
0.49 
5.11 
0.22 
0.22 
2.77 
2.23 
2.23 
8.31 



2.79 
2.78 
5.08 
3.07 
3.05 
3.41 
13.15 
12.35 
9.35 



2.59 
2.58 
5.23 
2.94 
2.92 
2.90 
13.10 
12.36 
8.80 



3.2 Asymmetries 

The forward-backward asymmetry (FBA) in the charged lepton decay is defined by 

ApB = Std(cos6 e )dr-S\d(cos6 e )dr 
Jq 1 d(cos8e)dT + J°_ 1 d(cos6e)dT 

where 0£ is the angle between the lepton £ + momentum in the dilepton cm. frame 
and the B s meson decay Lb. frame. In the SM, it is easy to show that the FBAs in 
B s — > (rj: f]')£ + £~ are equal to zero and the FBA in B s — > (f)£ + £~ is given by 

12s<4 (l - f )' (1 - r,) Re(G v F / + G° A Fy) 

In Figure 5, we display Afb(s) in B s — > <p£ + £~ (£ = fi, r) as a function of s. As seen from 
the figure, Afb(B s — > <fifi + fi~) changes sign at s = s ~ 0.1, i.e., Afb(s ) = 0, where s 
satisfies the relation 

s ° = Tw/" 2 ^ 1 + v^Htt + ( 35 ) 

Similar to B — > K*£ + £~ with £ = e or the location of s depends on the WCs and 
is insensitive to the hadronic uncertainties [5, 28]. 
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By writing the unit vector of n — ^ = ±1, we can define the longitudinal lepton 
polarization asymmetry in B s — > M s £ + £~ as follows 



dr(n=-l) 



dr(n=l) 



dr(n=-l) , dr(n=l) 
^/^ ~r 



ds (is 

From Eqs. (22), (23) and (36), P L in B s -> P^+r (P = 77,77') and P s 
by 



(36) 

4>£ + £~ are given 



= (c. 



P\2 



2(l"f) ! 



(l + f)«P + t5 P 



Re 



VP \ C e s ff F + - 2- 



C 7 F, 



7 r T 



pf = 



2(l - f I 

v S J R e (_ 

(l + f)a* + tfy l *> 



1 + y/rpJ 
G+F+* + ^(l-r )(l-r^- S ) 



(37) 



(G° Pf + G+P°*) + (^ + 12t> s)(l - r ) 2 G° F° A * + 8^ s GyPy] } (38) 

respectively, where C# = — sin 9 and Cg = cos 0. 

In Figure 6, we show the longitudinal lepton polarizations in B s — > (77, 77', <fi)£ + £~ with 
£ = H and r. The results for the electron modes are similar to the muon ones. We remark 
that, as shown in the figures, (s) are close to -1 for the muon modes except the end 
points of = 4m 2 and q^ax = ( m B s ~ m M) 2 in which they are zero and they range 
from —0.4 to and —0.2 to for the r ones, while the average ones of < P2 W) > without 
LD effects are — 0.97(— 0.96) and — 0.32(— 0.17) for the muon and tau modes, respectively. 
For B s — > <p£ + £~, P£(s min ) = but P£(s m ax) 7^ since cannot be factorized out in 
Eq. (38). We have that < P^ > M = -0.81 and < P^ > T = -0.49. 

4 Conclusions 

We have studied the rare B s decays of B s — > (77, 77', (j))vv and B s — > (77, 77', 4>)£ + £~ 
{£ = e, /x, r). In our analysis, we have used the form factors of B s — > (77, 77', 0) 
transitions calculated in the LFQM and CQM. We have found that Br(B s — > 
v ££) {£ = z/,e,/i,r)= (23.4,3.43,3.42,0.67) and (21.7,3.13,3.12,0.67) x 10~ 7 , Br(B s -> 
r]'££)=(25.2, 3.64, 3.63, 0.50) and (23.8,3.42,3.41,0.43) x 10~ 7 , and Br(B s -> 0£l) = 
(12.02,1.72,1.64,0.16) and (11.65,1.69,1.61,0.15) x 10~ 6 without LD contributions in 
the two models, respectively. We have also discussed the longitudinal lepton polarization 
and forward-backward asymmetries in the charged lepton modes. We have shown that 
the behaves of the asymmetries are similar to those in B — > K^£ + £~. Clearly, some of 
the above rare B s decays and asymmetries can be measured at the BTeV and LHC-B 
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experiments. 
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Figure 1: One- loop diagrams for the short-distance contribution in the standard model. 
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(a) (b) 




Figure 2: Form factors as functions of s — q 2 /M^ for (a) B s — > i] s (m v ) and (b) B s — > 
f] s (m' v ) with solid, dash, and dash-dot curves representing F + , F_, and F T , and (c,d,e) 
B s — > 0, respectively. The thick and thin lines stand for the results from the LFQM and 
CQM, respectively. 
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Figure 3: Differential decay branching ratios as a function of s = q 2 /m% s for (a) B, 
i] vu, (b) B s — > rj'uu and (c) B s — > <p vv. Legend is the same as Figure 2. 
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Figure 4: Differential decay branching ratios for B s — > i]£ + £~ in (a-c), B s — > i]'£ + £~ in (d- 
f), and B s — > <p£ + £~ in (g-i) with £=e, /i, and r, respecitvely. The curves with (without) 
resonant shapes represent including (non-including) LD contributions. Legend is the same 
as Figure 2. 
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Figure 5: Forward-backward asymmetries in the LFQM for (a) B s — > (f)fi + /j and (b) 
B s — > (f)T + r~ . Legend is the same as Figure 4. 



19 




Figure 6: Longitudinal lepton polarization asymmetries in the LFQM for B s — > r/£ + £~ in 
(a,b), B s — > r]'l + £~ in (c,d), and B s — > (f)£ + £~ in (e,f) with £=(/i, r). Legend is the same 
as Figure 5. 
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